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ABSTRACT 

We examine the properties of dark matter halos within a rich galaxy cluster using 
a high resolution simulation that captures the cosmological context of a cold dark 
matter universe. The mass and force resolution permit the resolution of 150 halos with 
circular velocities larger than 80 kms -1 within the cluster's virial radius of 2 Mpc. 
This enables an unprecedented study of the statistical properties of a large sample of 
dark matter halos evolving in a dense environment. The cumulative fraction of mass 
attached to these halos varies from 0% at 200 kpc, to 13% at the virial radius. Even at 
this resolution the overmerging problem persists; halos that pass within 200 kpc of the 
cluster center are tidally disrupted. Additional substructure is lost at earlier epochs 
within the massive progenitor halos. The median ratio of apocentric to pericentric 
radii is 6:1; the orbital distribution is close to isotropic, circular orbits are rare, radial 
orbits are common. The orbits of halos are unbiased with respect to both position 
within the cluster and with the orbits of the smooth dark matter background and 
no velocity bias is detected. The tidal radii of surviving halos are generally well-fit 
using the simple analytic prediction applied to their orbital pericenters. Halos within 
clusters have higher concentrations than those in the field. Within the cluster, halo 
density profiles can be modified by tidal forces and individual encounters with other 
halos that cause significant mass loss - "galaxy harassment" . Mergers between halos 
do not occur inside the clusters virial radius. 

Key words: cosmology: theory - dark matter - large-scale structure of the Universe 
- galaxies: clusters - galaxies: halos - methods: numerical 



1 INTRODUCTION 

Rich clusters of galaxies are large cosmological laboratories 
that may provide unbiased mixtures of the matter content 
of the Universe (White et al 1993, Evrard 1997). Clusters 
are prominent structures in the Universe; their evolution can 
be followed with samples out to z ~ 1 (Rosati et al 1998). 
They are the most massive virialised objects in the Universe 
and are the most recent objects to form in hierarchical for- 
mation models. Their masses can be determined by several 
independent methods (e.g. Carlberg et al 1997, Cen 1997, 
Wu & Fang 1997, Smail et al 1997). 

The evolution of the mass function of clusters is sen- 
sitive to key cosmological parameters (e.g. Bahcall, Fan & 
Cen 1997, Bartelmann et al 1997, Borgani et al 1997, Carl- 
berg et al 1996, Eke, Cole & Frenk 1997, Wilson, Cole & 
Frenk 1996). This evolution can be calculated by using ei- 



ther analytic methods (Press & Schechter 1974) or cosmo- 
logical N-body simulations (Eke et al 1997). The weakness 
of analytic methods is their inability to follow halos that 
accrete into larger systems. In the past, numerical simula- 
tions have shared this problem. To sample a large volume, 
the poor resolution within virialised systems leads to soft, 
diffuse halos that are rapidly dissolved by tidal forces. This 
is the classic over-merging problem (White et al 1987) that 
has lead to problems when comparing the mass distribution 
within dark matter simulations to the observed properties 
of galaxies (Summers et al 1995). 

We were determined to perform simulations that re- 
solved the distribution and evolution of galaxy halos within 
clusters. There are many questions that we will address with 
these simulations: What is the orbital distribution of the 
galaxies within clusters and are they biased in any way? 
What is the extent of galactic dark halos within clusters 
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and how much of the cluster's mass distribution is attached 
to galaxies? How do these properties evolve with time and 
within different cosmological models? Are the density pro- 
files of isolated "field" halos similar to the halos that form 
within the environment of a rich cluster? How does the clus- 
ter environment modify the internal structure of halos? The 
frequency of mergers between halos within the cluster en- 
vironment and the heating rate from halo-halo encounters 
are questions of importance for studies of the morphological 
evolution of clusters. These are amongst the many questions 
that have remained unanswered due to the over-merging 
problem. 

Some have suggested that the only way to avoid merg- 
ing is to follow the evolution of the baryons — even if they are 
only ~ 5% of the mass (Evrard et al 1994). This assertion 
seems counter-intuitive; it's now clear that mass and force 
resolution can overcome the overmerging problem in dark 
matter simulations (cf. Moore, Katz & Lake 1996, Klypin, 
Gottlober & Kravtsov 1997, Brainerd, Goldberg & Villum- 
sen 1997, Moore et al 1997). With fast parallel computers 
and highly tuned algorithms, our simulations have hundreds 
of surviving halos within the virial radius of a rich cluster 
(Moore et al 1997). Increased mass and force resolution lead 
to higher central densities in galactic halos, enabling greater 
survival within a cluster. For the first time we can compare 
the mass distribution with the galaxy distribution in a rich 
cluster. 

The plan of this paper is as follows: In §2, we describe 
the N-body simulation, techniques and parameters. In § ^ we 
describe two algorithms to identify "halos within halos" . Af- 
ter creating a catalogue of halos, § ^| turns to results on their 
global properties; sizes, masses, radial distribution, orbital 
properties and merger histories. § |B] discusses the internal 
properties of halos: their density profiles, correlations be- 
tween structure parameters and global parameters and the 
evolution of these quantities with time. We conclude in § El 



2 THE TV-BODY SIMULATION 

One of the goals of performing cosmological numerical sim- 
ulations is to compare the distribution and bulk properties 
of dark matter with the distribution and properties of the 
observable galaxies. A direct comparison has never been pos- 
sible, since structure in high density regions has been quickly 
erased as a consequence of numerical resolution - the over- 
merging problem. "Galaxies" are typically selected from the 
mass distribution in a cosmological simulation using a biased 
sample of dark matter particles. Previous studies of clus- 
ter substructure has been limited to using "galaxy tracers" 
(Carlberg 1991, Summers et al 1995), such as following the 
most bound particle of a halo before it becomes disrupted by 
the tidal field of the system. These hueristics enable the use 
of low resolution simulations, but their validity is anyone's 
guess at this point (Summers et al 1995). 

There are now several codes that are able to simulate 
a gaseous component. Although these codes are invaluable 
for many astrophysical problems; the original motivation be- 
hind these techniques was the hope of resolving galaxies in a 
cosmological context, thus solving the overmerging problem. 
They hoped to form galaxies and preserve them by increas- 
ing the central densities owing to gas dynamics. In turn, the 



halos would be more robust to disruption. Indeed, SPH sim- 
ulations of individual clusters, do give rise to a set of galaxy 
tracers that resemble a "real" cluster. However, the mass 
resolution in the dark matter component is not sufficient to 
resolve the dark halos of the galaxy tracers, typically one is 
left with a cold gas blob orbiting within the smooth clus- 
ter background (Frenk et al 1996). There are a variety of 
pathologies that arise if one uses too few particles to sim- 
ulate a large dynamic range in scales. Our simulations are 
designed to resolve scales of 5 kpc using dark matter. This 
would be a minimal resolution to simulate galaxy forma- 
tion with gas dynamics. One needs a fiducial dark matter 
simulation at high resolution to see differences owing to gas 
dynamics. One must insure that these differences are sensi- 
ble as an external check. With the gas representing <^10% 
of the mass, gross changes in numbers and orbits of galaxies 
would be surprising. 

Our aim is to achieve very high spatial and mass res- 
olution within a rich virialised cluster drawn from a "fair 
volume" of 100 Mpc 3 in a standard CDM Universe. In such 
a volume, there are several rich clusters, and none domi- 
nate the environment in an undesirable way. Any simulation 
method limits the number of particles that can be invested 
in a single simulation. Previous simulations of clusters suf- 
fered limitations due to the small volume used, forcing the 
run to stop at z ~ 1, (Evrard, Summers & Davis 1994), or or 
the use of vacuum boundary conditions outside the cluster 
(Carlberg 1991). 

With current technology, we can perform a single large 
simulation of ~ 10 8 particles or try to tackle a number of 
different problems using simulations with ~ 10 7 particles. If 
we simulate our "fair volume" at uniform resolution, there 
will be ~ 10 4 particles within the virial radius of a cluster, a 
resolution that is insufficient to resolve substructure. Previ- 
ous attempts to resolve the inner structures of cluster halos 
using ~ 10 5 particles failed to resolve more than a handful 
of satellite halos (Carlberg 1991, Carlberg & Dubinski 1991 
, Tormen 1997, Tormen et al). To achieve higher resolution 
within an individual cluster we initially perform a simulation 
of a large volume of a CDM universe as described above, nor- 
malised such that as = 0.7 and the shape parameter F = 0.5 
(H=50km s _1 Mpc -1 is adopted throughout). 

We used a nesting scheme that we call "volume renor- 
malization" to achieve higher resolution within a region of 
greater interest. This technique has been used to probe 
quasar formation at high redshift (Katz et al 1994) and to 
follow the density profiles of halos in a cosmological context 
(Navarro et al 1996). We generate initial conditions (ICs) for 
the volume at two resolutions, one that places ~ 10 7 parti- 
cles within the entire volume and one such that there would 
be <^10 6 particles in the targeted cluster. We run the lower 
resolution model and select a virialised cluster at z = 0. The 
particles within about twice the virial radius of the cluster in 
the final state are traced back to their locations in the ICs. 
Within this region, we use the higher resolution ICs. Beyond 
this high resolution region the mass resolution is decreased 
in a series of shells by combining particles in the high reso- 
lution ICs at their center of masses. In this way, the exter- 
nal tidal field is modelled correctly. The starting redshift in 
the high resolution run is increased to z=69 such that the 
perturbations in the smoothed density field of the high reso- 
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lution region obey the constraint 5p/p<0.1. We then re-run 
the simulation to the present epoch. 

We use a new high performance parallel treecode "PKD- 
GRAV" to evolve the particle distribution. PKDGRAV 
(Stadel et al , in preparation) has accurate periodic bound- 
aries and an open ended variable timestep criteria based 
upon the local acceleration (Quinn et al 1997). The code 
uses a spline softening length such that the force is com- 
pletely Newtonian at twice our quoted softening lengths. In 
terms of where the force is 50% of the Newtonian force, the 
equivalent Plummer softening length would be 0.67 times 
the spline softening length. In the high resolution region, 
our particle mass is 8.6 x 10 8 Mq. We perform two runs with 
10 kpc (RUN1) and 5 kpc (RUN2) softening lengths. The 
final virial radius of the cluster is ~ 2 Mpc and the mass is 
4.6 x 1O 14 M0 so that we have approximately 600,000 parti- 
cles within a sphere of overdensity 200. (Note that the clus- 
ter that is analysed here is the "Virgo" cluster from Moore 
et al 1997.) 
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3 CLUSTER SUBSTRUCTURE AND HALO 
IDENTIFICATION 

3.1 The density distribution in the cluster 

In past work, halos in dissipationless JV-body simulations 
have usually spontaneously dissolved when entering clusters. 
Two physical effects conspire with the finite numerical reso- 
lution to erase dark matter halos in clusters (Moore, Katz & 
Lake 1996). Halos are heated by cluster tides and halo-halo 
encounters, thus losing mass as they move into the poten- 
tial well. When the tidal radius approaches ~ 3 times their 
"core radii" (owing to either a density plateau or gravita- 
tional softening), they dissolve. Hence, it takes very high 
resolution to retain dark matter substructures at a distance 
100-200 kpc from the cluster's center. Our numerical pa- 
rameters were chosen so that halos would survive at these 
scales. The wealth of halos retained in our simulated cluster 
is visible in Figure jj] 

The upper panel is a map of the density distribution in 
a box of size R200 (see §3.2 for a precise definition), centered 
on the cluster and projected onto a plane. Each particle is 
colour coded according to the logarithm of the local den- 
sity (defined using an SPH smoothing kernel over 64 parti- 
cles in a code called SMOOTH [Stadel fc Quinn 19 97, http 
ref: http://www-hpcc.astro.washington.edu/toolsl). Only 



regions with density contrast S > 30 are shown. The clus- 
ter boundaries, set at -R200, correspond to the contours of 
the central blue region. Much of the mass inside ifeoo lies 
in the dark matter halos that we will analyze here. Their 
projected distribution is shown as a 'circle plot' in the lower 
panel of Figure |l|. The radius of each circle is the halo 'tidal 
radius' (Sec. 3.3) in units of J?200- Note that halos of simi- 
lar central densities (similar colour in the density map) may 
have largely different radii depending on their distances from 
the cluster's center. It is remarkable that substructure halos 
cover such a large fraction of the projected cluster's area. 



Figure 2. The density (upper panel) and circular velocity (lower 
panel) profiles for the cluster. The solid curves are at the final 
epoch and the dotted curves show the cluster at a redshift z = 0.5. 
The density is measured in units of the present critical density 
PO.cr', R is the physical distance to the cluster center, and was set 
at the location of the most bound particle. The arrow marks the 
value of the cluster virial radius -R200 at 2 = 0. 



3.2 Cluster properties and evolution 

We define the cluster center as the position of its most bound 
particle. This particle is within a softening length of the 
center of the most massive halo found by SKID. The density 
profile of the RUN1 cluster calculated in spherical shells is 
shown in the upper panel of Figure ^ (the solid line is for 
z = and dotted line for z = 0.5). The cluster forms at 
z ~ 0.8 from the mergers of many halos along a filamentary 
structure and, at z = 0.5, it has not yet virialised since it still 
has quite a lumpy structure but the global density profile is 
roughly similar to that measured 5 Gyrs later at z—0. (We 
shall compare properties of the substructure identified at 
both epochs.) 

The lower panel of this figure shows the circular ve- 
locity profile V C {R) = (GM(R)/R) 1/2 , where M{R) is the 
mass within R. The virial radius of the cluster is defined 
as the distance -R200 for which the average density enclosed, 
pc(R20o), is 200 times the cosmic density, p cr ; we obtain 
-R200 = 1.95 Mpc at z = and 1.2 Mpc at z = 0.5. The 
cluster is not spherical and has axial ratios that are roughly 
2:1:1. (In the following we will always use units of kpc and 
kms~ for lengths and velocities, unless we explicitly state 
otherwise.) 

Figure ^| shows the growth of the cluster's mass with 
redshift. Defining the formation redshift of the cluster as 
that where it has accreted half of its final mass, Zf orm ~ 
0.8 for our cluster. This is slightly earlier than expected for 
an average cluster of this mass from the Press-Schechter 
theory (1974; Lacey & Cole 1993), where Zf orm ~ 0.5 (as 
calculated from a routine kindly supplied by Paolo Tozzi). 
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Figure 1. Density map (colour plot) and circle plot of the halo radii within the cluster's virial radius (taken here as the length unit) at 

2 = 0. 
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Figure 3. Growth of the cluster's mass M c i with redshift z. M c [ 
is the mass encompassed by the cluster's 'viriaF radius at each z. 
The errorbars delimit the 1—cr confidence interval for the growth 
rate predicted using the Press-Schechter theory. 

This is not an unwelcome feature, since the cluster is in true 
virial equilibrium at z = 0. 

3.3 Halo identification 

Identifying DM halos in the high-density environment of the 
cluster is a critical step (cfr. Klypin, Gottlober & Kravtsov 
1997). The halos jump out visually, so while it is relatively 
straightforward to identify the halo centers, we must select 
only the bound particles to characterize the halo. We want 
to screen the cluster background that's streaming through, 
but the substructure itself will be tidally distorted and 
may have tidal tails of material that are loosely attached 
to the halo. Our group finding algorithm uses local den- 
sity maxima to find group centers and and then iteratively 
checks for self-boundedness to define group membership. 
Each group of particles found belongs to an individual halo. 
The algorithm is an improved version of DENMAX (ref), 
named SKI D and is fully described by Stadel et al (1996 , 
see http ref: jittp: / /www-hpcc. astro. Washington. edu/tooh ). 
For each simulation we adopt a linking length of 1.5l ao ft and 
a minimum number of member particles of 16, correspond- 
ing to a mass of ~ 1.35 • 1O 1O M0. In general, we shall use 
halos with more than 16 particles when they are employed 
as tracers, but we shall adopt a minimum number of 32 par- 
ticles when their individual properties are relevant. 

The high resolution region that we analyze is roughly 
the turn-around radius of the cluster or about twice the 
virial radius, Rt a — 2J?200 — 5 Mpc. Within this radius, 
there are 495 and 522 halos for RUN1 and RUN2, and 208 
and 227 halos respectively within i?2uo- Changing l so f t by 
a factor of two does not make much difference on a global 
scale, but if we restrict ourselves to the inner parts of the 
cluster the difference between the numbers of halos changes 
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Figure 4. Density profiles and circular velocity profiles are plot- 
ted for two large dark matter halos extracted from the simulation. 
The left panel is a "peripheral" halo that lies beyond the main 
cluster's virial radius at R = 4.2 Mpc. The halo in the right panel 
lies within the cluster at r = 1.2 Mpc. The radius r measures 
the distance from the center of each halo. In the latter case, the 
radius at which the high density background of cluster particles 
dominates the halo particle distribution is apparent as a flatten- 
ing of p(r) and the linear rise of v c (r) in the outer regions. We 
denote the peak value of the circular velocity as v pea k and the 
radius at which it occurs as r pea j.. 



significantly: Inside R < i?2uo/2, RUN2 has 91 halos com- 
pared to the 59 found in RUN1 and RUN2 has twice as 
many within R < ifeoo/4 (30 instead of 16). The innermost 
halos in RUN1 and RUN2 are at ~ 200 kpc and ~ 100 kpc 
respectively. These differences reflect the softening length's 
effect on the halos' central densities that determine their 
survival against tidal disruption (Moore et al 1997; see also 
the discussion related to Figure n below) . Quality control of 
our halo finding algorithm was insured by visually inspect- 
ing the density distribution inside ifeoo to verify that we 
neither missed nor created structures. 

We use the output of SKID to determine the halo struc- 
tural parameters. In particular, it estimates the extent of a 
halo using the distance to the least bound particle. How- 
ever, the full 6-dimensional phase space information is never 
available in the real Universe, therefore we shall compare 
results from SKID using an "observable" quantity for each 
halo. For example, the (projected) mass distributions can be 
determined using either weak lensing, the rotation curves of 
spirals or the velocity dispersion profiles of ellipticals. Moti- 
vated by gravitational lensing, one possibility is to define a 
halo's radius using its density profile p(r), where r is the dis- 
tance from a halo's center, and measuring the radius (= r p ) 
at which the local density of the cluster background domi- 
nates and p(r) flattens (cfr. fig. Q). Rotation curves or veloc- 
ity dispersion profiles of isolated objects will eventually de- 
cline with radius, but if a halo is embedded within a deeper 
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potential, at some radius its profile will turn around and in- 
crease as the velocity dispersion of the cluster background 
starts dominating. (Figure 4 shows that the position of these 
inflexion points are essentially equal.) We therefore combine 
these two definitions and use the inflexion point (= r vc ) of 
the effective circular velocity v c = (GM(r)/r) 1 '' 2 as an al- 
ternative independent estimate of the extent of a halo (cfr. 

The circular velocity is less noisy than the density or 
velocity dispersion and thereby more suitable for an auto- 
mated procedure. Moreover, the inflexion point of v c (r) can 
be easily detected by searching for a minimum, without any 
knowledge of the background density (as it would be nec- 
essary, if we were to implement an overdensity criterium). 
The radius r vc can either overestimate r p , for steeply declin- 
ing velocity profiles, or underestimate it, for profiles close to 
isothermal. If the halos are described by the NFW model 
(Navarro, Frenk & White 1996; see also Section 5), the for- 
mer condition applies to halos with v c ~ 50kms _1 and the 
latter to halos with v c ^, 200 km s -1 for background densi- 
ties ~ 300 times the cosmic average (as we will see however, 
the profiles of tidally 'stripped' halos decline with r more 
steeply than an NFW profile). In our case, for small halos, 
the difference can be at most ~ 10%; as for large halos, r vc 
differs from r p by i$ 20% in about 1/3 of the halos in our 
sample with 120 km s -1 . These differences do not sig- 
nificantly affect our results. With this definition, there is a 
contribution to the mass encompassed by a halo from the 
smooth background of order 20%, which we subtract from 
the quoted halo masses. 

For each SKID halo, we calculate p(r) and v c (r) using 
equally spaced bins of 2 kpc, such that the number of parti- 
cles in each bin is nearly equal. The departure of p(r) from 
isothermality is betrayed by a peak of v c (r), v p k occurring 
at r v k. The catalog values of r p k and v p k are estimated b; 
fitting a cubic spline to v c (r). The left panels of Figure 
show a sample p(r)/po, C r and v c (r) curves for a large pe- 
ripheral halo at a distance of R = 4.2 Mpc from the center 
of RUN1 (hereafter cluster halos are those within R200, and 
peripheral halos those between -R200 and Rto)- 

Halos beyond the cluster boundary can be easily char- 
acterized by their 'virial' radii and masses M200 = M(r2oo)- 
For isothermal spheres the circular velocity at V2oo is V200 = 
(r2oo/kpc)/ikms -1 to within a few percent. For the halo 
mentioned above, r 2 oo ^ 400 kpc (M200 = 3.4 ■ 1O 12 M ). 
The right panels of Figure ^Jshow p(r)/po, cr and v c (r) for a 
tidally limited massive cluster halo at R = 1.2 Mpc, where 
pikg — 300po,cr- The halo radius r na i is either the virial 
radius, V2oo or the tidaJ radius rud, as appropriate. Sim- 
ilarly, the halo mass is defined as Mhaio = M(rh a io) = 
Vc{rhaio) r haio /G. 

Force softening and the finite mass resolution introduce 
halo cores that are visible in Figure |as a flattening of p(r) 
at r£ 10 kpc. We define a core radius r core for each halo as 
the radius where v c has risen to 70% of v pea k- 

Once the halo positions and structural parameters are 
known, we can start to address the questions raised in the 
Introduction. We begin by examining the global properties 
of the halo distribution, with a preliminary discussion of 
numerical resolution effects. 
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Figure 5. For every halo within R = 5 Mpc we plot the ratio of 
halo core radii and outer radii vs. clustercentric position. Halos 
contain at least 32 particles. 



4 PROPERTIES OF THE HALO 
DISTRIBUTION 

4.1 Numerical resolution effects 

Halo in clusters are "harassed" by encounters with other 
halos combined with strong global tides. This adds internal 
energy and leads to mass loss. Using a static cluster poten- 
tial, Moore, Katz & Lake (1996) found that halos rapidly 
disrupt when r t idai < 3-4r core . Our fully self-consistent sim- 
ulations verify this. Figure ^ shows the ratios r core /r ha i vs 
R for our halo set using RUN1. This ratio increases towards 
the cluster center and no halos exist with r CO re/rhaio^ 0.3. 
(Similar results hold for RUN2.) Clearly, the softening will 
set a floor to the core radius, but in general we find that 
r CO re ~ ho ft, and correspondingly, that halos smaller than 
~ Shoft have all dissolved. The lower boundary for Thaio in 
RUN1 and RUN2 is indeed ~30 and ~15 kpc respectively. 

Numerical limitations can also be seen in Figure |(| a 
plot of v pea k versus R for RUN1. There is no bias apparent 
in this plot, halos of all circular velocity exist over a wide 
range in R. We show two curves that exclude regions of space 
owing to our numerical resolution: (i) defining a minimum 
particle number to identify halos sets a lower mass limit, 
Mi irn = 1.35 ■ 10 10 Mq, (ii) rtidai is correlated with R and 
halos dissolve when rudai^j 3r core ~ 3l so f t - Approximating 
the halos as isothermal spheres that are tidally stripped like 
layers off an onion, the mass and circular velocity are re- 
lated as: M ha io = M200 = Gv^r2oo for unstripped peripheral 
halos and Mhaio — Gv^rud- for "stripped" cluster halos. 
The tidal radius obtained from v c through rud — Rv c /V c 
(V c — lOOOkms -1 is the circular velocity for the cluster) in 
the approximation that a halo at R is tidally truncated "lo- 
cally" , 1. e. according to the value of pc at R (this is a limit- 
ing case as their pericentric radius can only be smaller). Ap- 
plying condition (i) leads to the dotted line in Figure p[ crite- 
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Figure 6. The peak circular velocities of halos, f pea fc, is plotted 
against their clustercentric position R at z = 0. The lines give the 
expected limiting R at which halos of given v peak can be resolved 
because of the finite spatial resolution (dashed curve) and mass 
resolution (dotted line). 

rion (it) leads to the dashed line (we set v c = const = v pea k)- 
If tides due to the cluster potential are the only cause 
of halo disruption, our sample should be complete for ha- 
los with VpeakfZj 80kms _1 and pericenters that have always 
been greater than Rum ~ 250 kpc. 

4.2 Spatial distribution of halos 

The number density of halos, rihaios, is plotted against R 
in Figure for the two runs. The dashed line refers to the 
particle density, n par t, normalised to the halo density in the 
interval [R200,Rta]- 

Beyond the cluster's virial radius the curves for halos 
and particles have similar slopes. Within the cluster's virial 
radius we see that the halo distribution is "antibiased" , i. e. 
less concentrated, with respect to the mass distribution. The 
halo number density profile is consistent, for R > 1 Mpc, 
with the average galaxy number density profile derived by 
Carlberg et al (1997) for the clusters of the CNOC survey. 
However, the average cluster mass profile that they derive 
from the same data is much shallower than that of the rel- 
atively small cluster analyzed here. 

If we regard the scaled particle distribution as the 
"asymptotic" halo distribution in the case of no bias and in- 
finite resolution (and sufficiently small physical cores), then 
we can estimate the numbers of missing halos by integrating 
the difference over R. For RUN2 we find that ~ 240 halos are 
"missing", about 50% of the halos with v pea k > 80kms _1 . 
The systematic difference between RUN1 and RUN2 makes 
tidal disruption of "softened" halos the most likely cause of 
this bias. However, the destruction of halos may not be due 
to numerical resolution alone. For example, binary merg- 
ers between halos of similar mass will lead to a single halo 



Figure 7. The number density of halos with v p< , a k > 80kms — 1 
as a function of R, for RUN1 (dotted line) and RUN2 (solid line); 
the former has a spline softening length of 10 kpc, twice that of 
the latter. The dashed line is the particle density scaled in such a 
way that, in the distance interval [i?200, Rta], the scaled number 
of particles equals the number of halos. The errorbars are 1— tj 
Poissonian errors from the counts in each bin; the arrow marks 
the value of i?200- 

with no memory of its history. This halo anti-bias may be 
overcome with higher resolution simulations that accurately 
resolve the structure of the smallest halos. 



4.3 Distribution of halo radii 

The extent of dark matter halos attached to galaxies in clus- 
ters has become directly observable via observations of grav- 
itational lensing (Geiger & Bartelmann 1997, Natarajan et 
al 1997). Here we can make some predictions for future sur- 
veys that will constrain the extent of halos as a function of 
position from the cluster center. 

The projected distribution of the halos within the clus- 
ter's boundaries was shown in the lower panel of Fig. |l| 
There is a clear decrease of halo sizes near the cluster's cen- 
ter even in this projected plot (there is little difference if we 
include halos up to twice ifeoo)- 

Figure ^ shows the average value of r halo I R200 as a 
function of clustercentric position, R/R200, at redshifts z=0 
and z=0.5. We use -R200 at each z as the length unit to 
highlight the self-similarity in the evolution of the cluster 
substructure. Halo radii clearly decrease as we move towards 
the cluster center, but the trend is hard to detect at 2 = 0.5 
because the cluster has accreted only relatively few halos (of 
Vpeak > 80kms _1 ), has a quite anisotropic mass distribution 
and tides have been efficient only in its very center. The 
mean size of halos beyond i?200 is ^ 8% of R200 and drops 
approximately linearly to zero as we move from the virial 
radius to the cluster center. (In Fig. ||, halo radii are those 
measured from the circular velocity profiles; the results using 
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Figure 8. In the upper panel, the average value of halo radii 
plotted against cluster-centric position at redshift z = (solid 
symbols) and z = 0.5 (dotted), for halos with i) pea fe > 80kms _1 . 
The lower panel shows the distribution against projected distance 
at z = 0. At each z, the unit is the cluster's virial radius at that 
epoch. The errorbars give the dispersion about the average. Note 
that the most massive cluster halo at z = has been excluded 
from the average; it has a radius of 470 kpc and is located at 
R = 0.9R 2 oo- 



Figure 9. The figure highlights the effects of tidal stripping on 
halo sizes at varying distance from the cluster's center (halos with 
v peak > 80kms — 1 ). For our halo sample, the ratio fhaio^/^peafc i 
with h = ho(l + z) 1,5 and units of kpc and kms" 1 , is approx- 
imately 70% of the value expected for a purely isothermal mass 
distribution. The dashed line gives the expected behaviour for 
isothermal halos being instantaneously stripped whilst falling into 
the deeper potential of the cluster. The upper and lower panels 
are for inwards and outwards moving halos respectively. See the 
text for further details. 



SKID radii are very similar.) As shown in the lower panel 
of this Figure, projection effects weaken but do not erase 
the trend of halo radii with R (here we have included halos 
up to 1.57?200, but it makes little difference changing the 
limiting R between once and twice ifeoo). Natarajan et al. 
(1997) do not detect a dependence on R, but presently, the 
observational uncertainties are large and the radial range of 
the data is limited. 

Halos of different v pea k obviously have intrinsically dif- 
ferent sizes. We can interpret Figure |^ and account for 
this intrinsic scatter using our "locally stripped isothermal 
onion approximation" (with v c = v pea k', see Section 4.1) 



and considering the ratios r\ = rhaio/{h~ 1 v pea k)- When ha- 
los form in isolation, r^aio = ?"200 ~ h~ 1 v c and r\ ~ const ~ 
1 ■ kpc/ kms -1 . Within the cluster, ruaio = rud — Rv c /V c ~ 
(h~ 1 v c )R/ 'Rwo) and r)(R) ~ R/R2oo- If halos have gone 
past pericenter r per i and their radii rud are determined by 
Pc (r peri), then rj w ill be smaller by a factor of r per i/R (cfr. 
Fig. in Sec. Q. 

In a fully virialized system there is no difference be- 
tween the motions away and towards the system's center. 
However, in a fi = 1 cosmology, clusters never stop accret- 
ing material and we might expect some different morpholo- 
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Figure 10. Same as in Figure ^ but the halo radii are measured 
using the SKID algorithm, that discards unbound particles, and 
data are shown for z = 0.5 too (lower panels). Only data for 
cluster halos are plotted; at z = 0.5, the formal cluster's virial 
radius is -R200 = 1200 kpc. 



gies between ingoing and outgoing halos. Moreover, halos 
reaching apocenters beyond the cluster's boundaries might 
be heated by low-speed encounters with other halos and re- 
expand, partly covering the effects of tides at the previous 
passage at pericenter. 

The distribution of 17 at z = is shown in Figure []. It 
has a trend similar to that expected for locally truncated 
isothermal spheres (dashed line). The separation between 
the latter and the points in the "periphery" (R > 2Mpc) 
is a measure of the departure of the actual density profiles 
of peripheral halos from isothermal. All outgoing halos must 
have necessarily passed their pericenters; therefore, 77 should 
depart from the dashed line more than those for infalling ha- 
los and we do observe evidence of this effect in the Figure. 
The corresponding results using SKID halo radii are shown 
in Figure hd which also shows results for z = 0.5 (only for 
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Figure 11. Fraction of mass within resolved halos in spherical 
shells of thickness 600 kpc versus cluster-centric position at 2 = 
0. The solid line is for RUN1 and the dotted line is for RUN2 
(slightly shifted for clarity). 



halos with R < R200 

zontal axis from z — to z 
-r 



note the different scale on the hori- 
0.5; h = ho(l + z) 1 ' 5 is H z 
in units of 100 Mpc/kms -1 ). The results for z — confirm 
the general picture illustrated by the previous Figure |J but 
here the points have a larger scatter and the trend with 
R is weaker, especially for outgoing halos: this is expected 
because, in the "dynamical" definition used by SKID, the 
tidal radii are less sensitive to pbk g (R)- At z = 0.5, there 
is evidence that tidal effects are already operating at this 
relatively high redshift, in agreement with the trend shown 
by Figure 



4.4 Fraction of mass in halos 

Figure [ll] shows the fraction of the mass attached to halos. 
As expected, it decreases sharply approaching the cluster 
center: it is % 5% at R^ 500 kpc and increases to ~ 20% 
at i?200- The total fraction of mass attached to halos within 
-R200 is about 13%. (These values do not depend sensitively 
on the adopted value of the softening parameter.) Outside 
the cluster, the halos account for about 20% of the total 
mass. The peak at R ~ 1.5 Mpc is not significant: it is caused 
by the largest halo within the cluster of mass 2.3 • 1O 13 M 
that contributes alone 5% of the total mass of the cluster 
and half of the mass in that bin. 

We can compare the mass fraction of peripheral halos 
with the Press-Schechter approximation (1974), using again 
a routine supplied by P. Tozzi. Adopting a minimum halo 
mass of 1.35-10 10 M© (16 particles) and the mass of the most 
massive peripheral halo (3.4 ■ 10 12 M Q ), we find ~ 0.25 from 
the analytical theory, in good agreement with the average 
of the last five outer bins, ~ 0.23. 
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Figure 12. Fractions of z=0.5 halos above resolution that have 
survived to z = retaining more than ~ 50% of their original 
masses; the errorbars are 1— c Poissonian uncertainties estimated 
from the numbers of halos in each distance interval. 

4.5 Merging and 'surviving' halos 

Comparing the distribution of the halos at z = 0.5 and 
z — 0, we can determine the merger rate and the fraction 
of halos that dissolve in this interval. We selected all the 
halos (of RUN1) with more than 32 particles at z = 0.5 
and examined their association with halos with more than 
16 particles at z = 0. We determine that a high-z halo 
is the progenitor of a low-z halo if the latter contains at 
least a significant fraction of the mass of the former, say 
25%. However, to account for mass loss via tidal stripping 
we need be less restrictive, therefore we considered also an 
extreme mass fraction cut-off of 1%, but our results are not 
very sensitive to this limit. 

We find, for halos with apocenters r apo > -R200 and 
at least one progenitor among those identified at z = 0.5, 
that 5-9% of them, depending on the mass fraction cut-off, 
are the product of a merger. There are a large number of 
mergers among halos with apocenters close to -R200 (3-5 out 
of 17). However, of the 38 halos with r apo ~ 80% -R200, none 
has merged. We find no evidence for merging once the halos 
have entered the cluster. This seems to argue against the 
possibility of Butcher-Oemler (1978, 1984) galaxies being 
cluster members with star formation switched on by mergers 
in the cluster environment (Couch et al 1994). 

Figure [12] shows the fraction of halos identified at 
z — 0.5 in different distance intervals that are also asso- 
ciated with z — halos. Although we adopt the extreme 1% 
mass fraction cut-off to establish this association, the results 
shown in the Figure essentially give the probability that ha- 
los will survive with 50% of their original masses attached. 
In fact, z — 0.5 halos with masses just above the limit of 
32 particles will be resolved at the 16-particle limit of the 
z — sample, only if they have not lost more than 50% of 
their masses; and these small halos dominate the statistics. 



At 7?iS 0.5 Mpc, such probability is ~ 0.5, against a value 
of order unity in the cluster periphery. 

4.6 Orbital parameters 

We now shift our attention to the motions of the halos within 
the cluster. Halos that follow radial orbits are more likely 
to be disrupted than halos that have circular motions, since 
the former penetrate further into the cluster potential well. 
Do we detect a bias in favour of circular orbits? We calculate 
orbits by approximating the cluster potential as a spherical 
static potential, computed using the density profile at 2 = 
and then using the position and velocity information for 
the halos. For comparison, we also compute the orbits of a 
random subset of 20,000 particles. 

In the upper panel of Figure ^ we plot the ratios 
r P eri/r a po vs R for the halos of RUN2 (RUN1 gives iden- 
tical results). The lower panel shows the average values of 
this ratio for the halos (solid line) and the particle sample 
(dotted line). The errorbars show 1-a standard deviations. 
We find that radial orbits are quite common and circular 
orbits are rare. The median ratio of apocentenpericenter is 
approximately 6:1 and nearly 25% of the halos are on orbits 
more radial than 10:1. A rough calculation by the authors 
reveals that this is very close to an isotropic orbital distri- 
bution within an infinite isothermal potential. 

We do not detect a large difference between halo or- 
bits as a function of ii, nor do we find a difference between 
the orbits of the particle background and the halos. This 
is surprising since we expected to find fewer halos on ra- 
dial orbits near the cluster center. The expected bias could 
be disguised if the central overmerging problem originated 
within the dense clumps that formed before the final cluster. 
Finally, we note that the radial velocity dispersion of the ha- 
los within the cluster is 720 km s -1 , a value that is within a 
few percent of the dispersion of the background particles - 
even when only the most massive halos are considered. 

Information on the distribution of the orbital param- 
eters of halos can be used to model the effects of tidal 
stripping and dynamical friction for halos within halos in 
semi-analytic models of structure formation based on the 
Press-Schechter approximation. Figure and |l| plot the 
probability density distributions of "circularities" and peri- 
centers for cluster (solid line) and peripheral (dashed) halos. 
For each halo, the "circularity" e = J/Jc(E) is defined as 
the ratio of its angular momentum to that of a circular or- 
bit with the same energy (Lacey & Cole 1993). There are 
no marked differences between the two groups of halos, al- 
though the orbits of cluster halos are more close to circular 
and penetrate further into the cluster than those of periph- 
eral halos. Among the latters, 15% have pericenters outside 
the cluster's boundaries and 9% come as close as 200 kpc 
(O.I-R200) to the cluster's center. This condition is twice more 
frequent among cluster halos; in the whole sample the frac- 
tion of r pEr i < 200 kpc is 14%. These results are in good 
agreement with results presented by Tormen (1997). 

For such radial orbits, we expect that the tidal radii of 
the halos are determined primarily by the global tidal field 
of the cluster. We can check if this is correct by estimating 
the truncation radius at the pericenter of each halo using 
r tid ~ r perl ■ v peak /V c (the use of V C (R) at R = r peri , in- 
stead of the constant value V c = V200 = lOOOkms -1 does 



© 1998 RAS, MNRAS 0, 1-7 



Dark Matter Halos within Clusters 11 



1 

0.8 
0.6 
0.4 
0.2 




z = 



o o 



O o 
o 



9,5 



' ^ Q 




° o o 8 - 



0.6 



0.4 



0.2 



2000 4000 
R/kpc 

z = 



halos 
particles 



! I I l_ 

2000 
R/kpc 



4000 



Figure 13. The upper panel shows the ratio pericentric and apoc- 
entric distances for each halo plotted against its current position. 
The lower panel shows the average value of these ratios for the ha- 
los (dashed curve) and a random subsample of cluster dark matter 
particles (solid curve). The errobars are the dispersion about the 
average. 



not make any detectable difference, because the variation of 
V C (R) is £ 10% in the range 0.05 < R/R200 < 1 ). We can 
test this prediction for our outgoing halos that must have 
passed pericenter recently, enhancing the likely validity of 
our approximation. In Figure |l6| we plot the expected tidal 
radius, according to the above formula, against the value 
measured, for both methods (from SKID and from the cir- 
cular velocity profile v c (r)). The agreement is excellent for 
the SKID values, with the exception of the points marked as 
stars. These latter points represent halos that are on very ec- 
centric orbits such that r per i is less than 300 kpc: These halos 
are more likely to suffer impulsive collisions with other halos 
as they pass close to the cluster center. We note that many 
of these halos have tidal tails that may cause the measured 
tidal radii to be overestimated. The scatter in the correlation 
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Figure 14. Probability density distribution of halo "circulari- 
ties" e = J/Jq(E) (see text) for cluster (R < R200; solid line) 
and peripheral halos (dashed). 
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Figure 15. Probability density distribution of halo pericenters 
r P eri (line styles as in Fig. fll| ). 



increases when we measure halo sizes using the halo circular 
velocity profiles, but the trend is still apparent. The radii 
estimated with this method are more sensitive to tidal tails 
than SKID radii, since the latter can discard the unbound 
streams of particles. This is the reason why some points in 
the lower panel of the Figure correspond to larger rud than 
those in the upper one. 

Recalling Figures ^| and [n], we can now check if tidal 
stripping is indeed responsible for the low values of r\ = 
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Figure 16. We plot the measured values of halo tidal radii against 
their expected values, assuming that the halos have isothermal 
mass distributions that are tidally stripped at their pericentric 
positions. The upper panel uses halo tidal radii measured using 
SKID, therefore only self-bound particles are included to deter- 
mine halo sizes; the 'v c ' method uses halo circular velocity profiles 
and includes the cluster background. The points represent outgo- 
ing halos (at z = 0) with R < O.8-R200; the stars denote those 
with r per i < 300 kpc. 



Thaioj (h~ 1 v pea i e ). Using the information on r per i, we can 
'correct' the values of 77 for outgoing halos as 77 — » 77- R/r per i. 
The effect of this correction is shown in Figure [l^, in which 
the open circles are the original points of the upper-right 
panel of Figure [lo] and the filled circles are 'corrected' As 
expected, the corrected 77 scale as R/R.200- In conclusion, 
the isothermal model predictions for the tidal radii work 
well when the pericentric positions of the halos are used. In 
turn, this good agreement confirms that our estimates of the 
orbital parameters for cluster halos are correct. 
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Figure 17. The dashed line gives the expected dependence on 
R of the ratios r\ = T^^/vp,,^ for isothermal halos that are 
tidally limited by the cluster at a distance R as they fall into the 
cluster for the first time. Such ratios attain a minimum when halos 
reach pericenter and are maximally stripped; therefore, moving 
away from the cluster's center, outgoing halos have values of r\ 
systematically below the line (the open circles; same as in the 
upper— left panel of Figure If we 'correct' their 77s using the 
information on their pericenters, the points are brought back to 
the line (the filled circles). 



5 HALO INTERNAL PROPERTIES 

We shall now examine the internal structure of halos in more 
detail with the main objective of studying how the cluster 
environment affects them. In particular, we shall focus on 
the distribution of r pea k and v pea k, that provide informa- 
tion on their internal concentrations. Then, we shall study 
in detail a small subsample of well resolved halos and ex- 
amine their density profiles and the compatibility with the 
analytic fits generally adopted in the literature to describe 
dark matter halos (Hernquist 1990; Navarro, Frenk & White 
1996). 



5.1 Statistical distribution of halo properties 

A detailed analysis of isolated cold dark matter halos in iV- 
body simulations has been carried out by Navarro, Frenk & 
White (1996 and 1997, NFW in the following). Their sim- 
ulations had a mass resolution such that individual halos 
contained of order 10,000 particles and force softening that 
was 1% of the final virial radii. Over a wide range of masses, 
NFW found relations between the properties of their pro- 
files, r pea k and v pea k and their virial masses M v . Further- 
more, the density profiles of halos could be fit by a univer- 
sal formula with varying concentrations, c, that can be pre- 
dicted from their masses within a given cosmological model. 

Here we want to address two questions: 

(i) Do the same relations found by NFW for isolated ha- 
los hold for the peripheral halos that surround the virialized 
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cluster? These are relatively isolated, but they evolved in an 
environment that is special owing to the nearby cluster. For 
instance, the streaming motions it induces could have de- 
termined peculiar merging histories of the peripheral halos, 
and merging may play a crucial role in shaping their internal 
structures (cfr. Syer & White 1997; Moore et al 1997). 

(ii) Do the same relations apply for cluster halos, that 
are affected by tidal stripping and halo-halo encounters? 
Cluster halos could also have peculiar properties reflecting 
the high-density environment in which they formed. 

These two factors could then affect the properties of 
our halo sample in such a way that, for instance, the ratios 
v p k/r p k, that reflect the concentrations of the halos, are sys- 
tematically higher than for the field population (e.g. cluster 
halos form earlier or are biased towards high concentrations 
such that they have survived tidal disruption) Alternatively, 
impulsive mass loss may cause halos to subsequently "relax" 
and re-expand towards an 'equilibrium' configuration. This 
could be detected as a change in concentration versus orbital 
position. 

We shall focus on the following two relations, for clus- 
ter and peripheral halos, M ha io vs. v peak and r peak vs. v peak . 
The former has the interest of relating masses to velocities, 
that are in principle more easily observable, whilst the lat- 
ter gives information on halo concentrations and the related 
issues discussed above. 

The distribution of halo masses vs. v peak is illustrated 
in Figure [l^, for peripheral (upper) and cluster halos (lower 
panel) and compared with the relation found by NFW 
within the cosmological model we have adopted (NFW 
1996; dashed line). (For clarity, only points with v peak > 
100 km s -1 are shown, since the noise is large for small ha- 
los). In each panel, the stars mark inner halos of each group 
(see the caption of the Figure for details; note that some in- 
ner halos within the present "periphery" sample have highly 
radial orbits and were inside the cluster at an earlier epoch.) 
For peripheral halos, the agreement between our data and 
the NFW curve is excellent (note that most of the scatter 
owes to points with v pea k < 150 km s -1 that may be affected 
by the numerical noise). Cluster halos have a much larger 
scatter and the fit is shallower. Their masses are smaller for 
a given v pea k, as naturally expected for stripped halos and 
this is most noticeable for the inner cluster halos. 

The distribution of v pea k and r pea k is shown in Fig- 
ure again for peripheral halos in the upper panel and 
cluster halos in the lower panel. The dashed line gives the 
relation for the halos studied by NFW and the dotted line 
in the upper panel is a power-law fit to the halos with 
v pe ak > 150 km s -1 (that appear to be less affected by noise) . 
The behaviour of the peripheral halos clearly agrees well, 
within the scatter, with that expected from the analysis of 
NFW. On the contrary, cluster halos lie significantly below 
the dashed line: they are skewed towards smaller values of 
r pea k for a given v pea k , with a more prominent effect for the 
innermost halos. 

At z — 0.5, we observe a behaviour similar to that at 
z — 0, as is shown in Figures and ^l| (in these figures, 
the dotted line is again a fit to the points with v pea k > 
150 kms -1 and the dashed line, the NFW's relation at z = 0, 
is drawn for comparison). 

In conclusion, the high-density regions of the cluster 
and the "accelerated collapse" of the cluster halos can affect 
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Figure 18. Distribution of halo masses versus v pea k for 'periph- 
eral' (upper panel) and cluster halos (lower panel) at z = 0. The 
stars mark inner points, with R < I.5R200 an< I R < O.6-R20O) re- 
spectively for the two groups. The filled circle gives the values for 
the entire cluster. The dashed line is the power-law fit obtained 
by Navarro et al (1996) for isolated halos. The dashed lines are 
similar power-law fits to the data of the form M = 10@v a , where 
8 is the rms scatter in mass about the fit. 



not only the sizes of the halos but also their internal struc- 
tures. We examine this issue further in the next section by 
looking in detail at the properties of a sample of large well 
resolved halos. 



5.2 The large halo sample 

We now restrict our analysis to several well resolved, massive 
halos that are not as affected by force softening or mass res- 
olution, by following the evolution of those halos that have 
r P eak > 25 kpc and v peak > 150 km s -1 (from Figure |l9"[ 
we can see that these values delimit a region in the plane 
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Figure 19. Distribution of v pea ^ and r pe(1 j, for peripheral (upper 
panel) and cluster halos (lower panel) at z = 0. The points have 
the same meanings as in Figure |is|. As before, the dashed line 
gives the expected relation for isolated halos (Navarro et al 1996) 
and the dotted line is a fit to the points with v pea k > 150 km s — 1 . 



(v pS ak,r psa k) where the properties of halos behave quite reg- 
ularly) . 

To highlight differences between "cluster halos" and 
"peripheral halos", we selected two groups in the distance 
ranges R/R200 < 2/3 and 4/3 < R/R200 < 7/3. Again, 
these values of R correspond to those roughly separating 
different halo behaviours in Figure [l9| These selection crite- 
ria yield 7 halos in the first distance range (although names 
like "Yogi", etc... would be preferable, for brevity's sake, 
we shall simply call them: A, B, C, D, E, F) and 7 in the 
second one (G, H, I, L, M, N, O). The first 3 of each group 
lie in the inner part of their distance ranges (f?i$ 0.8 Mpc 
and 2.8 Mpc ^ R ^ 3.5 Mpc, respectively). All of these halos 
contain at least 1000 particles within rhalo- We excluded halo 
"O" because it is the product of a recent merger (z ~ 0.5). 



Figure 20. Symbols and lines are the same as in Figure |l8| but 
for z = 0.5. Here, R200 is the formal cluster's virial radius at that 
epoch. 

All of the other halos have only one progenitor at z — 0.5 
and the inner halos are well defined up to z ~ 2. 



5.2.1 Evolution of the profiles 

We now examine the evolution of the density profiles for 
this sample of halos. We consider four values of the redshift: 
z = 0, z = 0.5, z = 1, and z = 1.8. In our sample, all 
peripheral and outer cluster halos show evidence of major 
mergers at z ~ 1, but the inner cluster halos (A to D) have a 
well defined major progenitor up to z — 1.8 (although they 
have captured many small satellites and experienced minor 
mergers since z ~ 1). For peripheral and outer cluster halos 
we limit our analysis to z — 0.5, while for the innermost 
ones we will examine the profile data up to 2 = 1.8. 

Figure |22| shows the evolution of p(r) for the cluster 
halos. All the curves for the present epoch flatten at large r 
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Figure 21. Symbols and lines are the same as in Figure [19] but 
for z = 0.5. Here R200 = 1.2 Mpc is the formal cluster's virial 
radius at that epoch. 

where the smooth particle background density of the cluster, 
pbkg, starts dominating (halos from left to right in each row 
of figures have increasing distances from the cluster center). 
Such flattening also appears in the profiles of the z — 0.5 
progenitors of A, B and C because they are within the high 
density environment of the forming cluster at that epoch. 

The features of p(r) are highlighted in the correspond- 
ing circular velocity profiles v c (r) shown in Figure |%| Note 
that the location and height of the peak of v c (r) have both 
changed for the inner halos (upper panels) , whilst the pro- 
files of the outer halos E and F, in the lower row, are re- 
markably stable, apart from the change in pbkg- 

Halos A and D show the most significant evolution and 
have quite steep outer slopes. Halo A formed at the inter- 
section of the two filaments whose collapse originated the 
cluster and its structure has been heavily disrupted by the 
tidal field there, causing it to lose a huge amount of mass. 
Halo D formed in the outer parts of one of the filaments, and 
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Figure 22. Evolution of the density profiles of a sample of mas- 
sive cluster halos after their last major mergers. The solid, dotted, 
short-dashed and long-dashed lines are the profiles measured at 
2 = 0, 0.5, 1, 1.8 respectively (some halos have their last major 
mergers at z > 1 and correspondingly, not all lines are shown; 
the 2 = 0.5 profile of D could not be measured because this halo 
was too close to the cluster's center at that epoch). The radius 
r is the distance from the halo centers in physical kpc and p is 
measured in units of today's critical density po.cr- 
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Figure 23. Evolution of the circular velocity profiles for the 
same halos as in Figure |22| plotted using the same line types at 
the different redshifts. The velocities are measured in physical 
kms -1 . 



© 1998 RAS, MNRAS 0, 1-7 



16 Ghigna et al 



100 200 300 100 200 300 100 200 300 



I 1 1 1 L-f-l i ■ i i | ■ ■ i i | ■ i i 1 1 ■ 




1 1 I 1 1 1 I 1 1 I I I 




100 200 300 100 200 300 100 200 300 

r/kpc 



Cfi 

O 



1.4 



III 

- 




i \ y 1 1 


- 




/ 

/ A 






/ 

/ 

/ 


- 

— 


/ 

* /7 


n * 

U i 






•' 


F 


E 


/ 


/: 

• / / 

: / : 


• 

* * 




'{ • 

- '' \ B 




# ""a 


/ » 






/ 

III 


•c 


1 1 1 1 1 


2.2 


2.4 


2.6 2.8 


Lo; 


= peak' 


/ km s _1 



Figure 24. Evolution of the circular velocity profiles of a sample 
of peripheral halos after their last major mergers. The dotted and 
short-dashed lines correspond to the redshifts z = 0,0.5,1. (At 
2 = 0, halo G has a distance from the cluster's center ~ 4/3i?200 
but was inside the cluster at earlier epochs. 



Figure 25. Evolution of the cluster halos of Figure |22| in the 
plane (v pea i 1 ,r pea i c ). The positions at z = are plotted as filled 
circles, those at z = 0.5 as stars (not for D) and at z = 1 as trian- 
gles (only for A-D). The dashed line is the relation expected for 
isolated halos at z = from Navarro et al (1996) as in Figure Hjjl 



its fate has been similar. Halos B and C have also evolved, 
but the changes have not been so dramatic. The values of 
Vpeak have decreased however by ~ 15% between z = 0.5 
and 2 = 0. The remaining halos E and F, that do not 
show signs of significant evolution, formed in the periph- 
ery, but at z = they are moving outwards. However the 
pericenter of F is close to its present distance (~ 1.3 Mpc). 
That for E is ~ 600 kpc, yielding an expected tidal radius of 
~ 150 kpc; this agrees very well with the SKID value, but is 
25% smaller than the tidal radius measured from the den- 
sity profile (directly from p(r), or from v c (r)). This is due 
to the anisotropic distribution of mass around E, that ap- 
pears stretched along the orbit with a "trail" of particles on 
both sides. This tidal debris were contained within its virial 
radius at z = 0.5 and, though stripped during the passage 
at pericenter and presently unbound to it, are still moving 
closely apparently on the same orbit. The lack of evolution 
in the density profile of E is then due to the assumption 
of spherical symmetry and the contribution from unbound 
particles nearby. 

Finally, Figure ^ shows the corresponding v c profiles 
for the 6 peripheral halos considered. Except for case G, 
these halos have never entered the cluster and their profiles 
show no significant differences between z = and 0.5. Halo 
G formed in the periphery but has been orbiting through 
the cluster since z ~ 0.5 and has clearly lost mass. 

In summary, of the large cluster halos identified at 
z — 0, only those with -R200/2 have density profiles 
significantly different from those they had before the forma- 
tion of the cluster. However, one of the large halos in the 
outskirts of the cluster has survived a passage (or more) at 
pericenter and has evolved considerably too. In particular, 



the values of v pea k can change: ~ 10% in three cases, but 
20-25% in two other cases. 

It is interesting to examine how the halos "move" in the 
plane (r pea k, v pea k) of Figure jlil so as to determine if their 
evolution (under tidal stripping or halo-halo encounters) is 
the reason for the systematically lower r p k for the cluster ha- 
los with respect to field halos of same v p k. Figure |25| shows 
the evolution of r p t vs v p t for the 6 cluster halos consid- 
ered previously from z = 1 (when available) through 0.5 
to z — (from triangles to stars to circles in the figure). 
From these data, it seems that the points did not move sig- 
nificantly away from the NFW's curve. They move parallel 
to the line rather than downwards and those furthest from 
the line are also those that formed the nearest to primor- 
dial high-density regions (A and C, that formed near the 
cluster's center, but also D, that formed in one of the giant 
filaments that merged into it). This could be evidence that 
halos form earlier in high-density regions and are thus more 
concentrated than those forming in the field. This impres- 
sion is strengthened by Figure bq, that shows the evolution 
of the "concentration" v pea k / r psa k (the symbols are as in 
Figure ^). The points on the left are for the halos form- 
ing near the center (A and C), those on the right for those 
forming in the periphery (B, E, F); the intermediate case D 
is duely in the middle. Although the ratios change, the two 
groups are clearly separated at both epochs with higher val- 
ues of concentration for the former and lower values for the 
latter group. Recently, Lemson & Kauffmann (1997) have 
studied correlations between halo properties and local en- 
vironment using large A^-body simulations. Except for the 
mass distribution of halos, they do not find any such cor- 
relation. It should be noted however that their definition 
of locality is based on a scale of 5-10ft -1 Mpc, that is still 
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Figure 26. Evolution of "concentrations" «p ea fc/ r pea*; f° r the 
halos of Figure ^ and the points represent the same redshifts 
as before. The halos are ordered from left to right according to 
the overdensity of the region where they formed: A and B formed 
at the intersection of the filaments whose collapse originated the 
cluster, D along one of the filaments and the others in the "field" . 



Figure 27. Comparison between the density profiles (at 2 = 0) 
of a sample of peripheral halos (solid lines) and the NFW model: 
the dotted lines are the expected NFW profiles (from halo virial 
masses) and the dashed lines are fits of the NFW profile to the 
data. 



mildly non-linear today. The correlation we find seems to 
involve scales of a few Mpc that become non-linear at early 
epochs. It will be interesting to examine the significance of 
our result with larger samples of cluster halos within simu- 
lations of even higher resolution. 

5.2.2 Analytic fits 

The evolved density profiles of isolated halos in TV-body 
simulations are well described for a large range of masses 
by the analytic model of NFW (although increasing the nu- 
merical resolution causes steeper inner profiles, Moore et al 
(1997); this is not an issue here since our sample of halos 
have similar resolution as those in the NFW simulations): 

P(r) = 5c 

Per (cr / r 200) (1 + cr/r 200) 2 



(5-1) 



with S c 



200 



(5.2) 



3 (Zn(l + c) - c/(l + c)) ' 

Navarro, Frenk & White (1997) developed an analytic proce- 
dure that gives c as a function of the halo mass M200 in any 
hierarchical cosmological model, based on the expected red- 
shift of collapse of a density perturbation of mass M in the 
Press-Schechter (Press & Schechter 1974) formalism Here 
we address the following two questions: (i) does the NFW 
profile and predicted c provide a good description of the pro- 
files of our peripheral halos? As mentioned previously, this 
question is not trivial since the environment within which 
these halos evolved is perturbed by the intense gravitational 
field of the cluster; (ii) How much do the tidally stripped 
cluster halos depart from the NFW predictions ? For exam- 
ple, a steeper outer slope may be typical of these halos. For 



this reason, we also consider the Hernquist profile (Hernquist 
1990; HER in the following): 



Per 



(cr/r 2 oo)(l + cr/r 2 oo) 



(5.3) 
(5.4) 



with 5 C = —c(l + cy . 

This profile has the same inner form as NFW, but asymp- 
totes to r~ 4 on large scales instead of r~ 3 . 

First we consider the sample of large peripheral ha- 
los described above. Figure |2t] shows their density profiles 
p(r)/p cr fi- In each panel, the dotted and dashed lines are 
NFW profiles predicted using the measured M200 according 
to the NFW procedure mentioned above, and tit to the data 
using a standard \ 2 minimization technique with c as free 
parameter (using r2oo measured from the data); for the fits, 
we used the points in the radial range delimited by l so ft and 
?"200- (Two-parameter fits, with both c and r2oo as free pa- 
rameters, yield very similar results; we use two-parameter 
fits for the cluster halos below, since r 2 oo cannot be de- 
fined for them.) We plot the residuals between the data and 
the NFW profiles and Hernquist profiles in Figure Each 
panel shows the fractional difference (p mo dei(r) — p(r))/p(r) 
as a function of r, where p(r) are the data and p m odei(r) 
is the density corresponding to one of the analytic profiles 
considered. With the exception of halo G, the residuals for 
the NFW profiles are within ~ 20-30% at distances from 
the halo centers 15kpc At smaller distances, the profiles 
of our halos are steeper than the NFW curves with residu- 
als typically in excess of 30% for the fits and 40% for the 
predicted curves. Note that, for the same halo, the different 
NFW curves have concentrations that can differ by vari- 
ous amounts, typically from c ~ 13-14 (expected) to 16-18 
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Figure 28. Residuals (om war (r) — o(r)) I o(r) between the NFW 
profiles and density profiles of Figure t27| (dotted and dashed lines 
are again for the expected and fitted NFW profiles). The residuals 
of fits of the Hernquist profile to the data are shown as thin solid 
lines. The arrows mark the halo virial radii. (For clarity, we show 
1 data point every 5). 



Figure 29. Comparison between the density profiles (at 2 = 0) 
of a sample of cluster halos (solid lines) and fits of NFW (dot- 
ted lines) and Hernquist's (dashed lines) profiles. An arbitrary 
isothermal profile is drawn for comparison (thin solid line). 
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(fits). The data do not discriminate significantly between 
NFW and Hernquist's profiles. Only for halo G, as expected 
from its steep profile, the HER fit fares better, although it 
still has positive residuals of ~ 40-50%. 

Let us now examine the cluster halos. Their density 
profiles (again solid lines) are shown in Figure together 
with NFW (dotted) and HER (dashed) fits (for the latters, 
we use the data in the radial range delimited by l so f t and the 
value where p(r) flattens approaching pbkg)- Over the scales 
of interest, the two fits do not differ significantly and both 
underestimate the central concentrations of the halos. As 
before, this can be seen better by examining the residuals in 
Figure At intermediate distances from the halo centers 
(r^S 40kpc), both fits have negative residuals in excess of 
~ 20%, and in excess of ~ 30% on smaller scales. At larger 
distances, up to about 75% of the halo tidal radii, they fare 
generally well, with residuals well below 20%. The exception 
is halo A, that has, as halo G, a particularly steep profile: for 
this halo, HER is still an acceptable fit, while NFW has large 
residuals, ~ 40-50%. When r approaches rud, the fits depart 
rapidly from the flattening p(r), with average residuals at 
rud of order of 30% or more, for NFW, and about 10% 
higher for HER. 

As a general remark to conclude this section, if we limit 
the halo sample to lie within r < 3/4rud, Hernquist's fits 
are better than NFW for the profiles of cluster halos (or 
stripped halos in general), while the latter model fares well 
for peripheral halos. Although our resolution is not good 
enough to address reliably the issue of halo concentrations, it 
is remarkable that our halos are systematically more concen- 
trated than the NFW model predicts, when from the poorer 
resolution we would expect just the opposite. In our anal- 
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Figure 30. Residuals (pFlT( r ) ~ p( r ))/p( r ) f° r the fits of Fig- 
ure |29] (NFW, dotted line, and Hernquist, dashed). The arrows 
mark the halo tidal radii. 



ysis, the NFW profiles fit to the data have concentrations 
about 25% higher than expected and still, they generally 
underestimate the profiles measured at small scales. 
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6 SUMMARY AND DISCUSSION 

We have explored the consequences of increasing the force 
and mass resolution within a dark matter only simulation of 
a rich galaxy cluster that forms hierarchically within a cold 
dark matter simulation of a closed Universe. By resimulat- 
ing regions of interest using increased resolution, we have at- 
tained an unprecedented view of the internal structure of a 
massive dark halo. With approximately one million particles 
within the virial radius and force softening that is 0.25% of 
R200 , we resolve 150 halos with circular velocity larger than 
80 kms^ 1 within -R200 at z = 0. Most of these halos have 
made several orbits within the cluster and are easily identi- 
fied as potential minima or density enhancements above the 
background. 

This work demonstrates that the overmerging of dark 
matter substructure within virialised structures can be 
greatly reduced given high enough numerical resolution. The 
statement that "at bottom the problem [of overmerging] ap- 
pears to be not numerical but physical: gravitational dynam- 
ics alone cannot explain the existence of galaxy groups and 
clusters" (Weinberg, Katz and Hernquist 1997) is, at the 
light of these new results , completely wrong. 

Overmerging within the cluster environment is due to 
the disruption of halos by the global tidal field and halo- 
halo encounters (c.f. Moore et al 1996), probably primar- 
ily within the large dense halos prior to the formation of 
the main cluster. Although our softening length (5 kpc) is 
a small fraction of the cluster's virial radius, the rotation 
curves of the halos peak at radii of ~ 30 kpc. Hence the 
cores of these halos can be softened to a degree that ef- 
fects their evolution. This is the primary reason the halos 
are still being disrupted; given high enough force and mass 
resolution, it should be possible to overcome most of the 
overmerging problem within CDM simulations. 

The aim of this paper has been to analyse the properties 
and dynamics of the dark matter substructure and we find 
the following key results: 



(less centrally concentrated) with respect to the global mass 
distribution. 

• The density profiles of a sample of well resolved halos 
indicate that those forming in the high-density regions of 
the collapsing cluster have higher concentrations than those 
found in isolated environments. We show that this is most 
probably due to their earlier collapse redshifts rather than 
the internal response of the halos to mass loss and heating 
from tidal stripping. 

• Most of the halos within the cluster and in the clus- 
ter proximity have density profiles that are well fit by NFW 
profiles (Navarro, Frenk & White 1996). Halos that lose a 
great deal of mass through tidal stripping have outer den- 
sity profiles as steep as p(r) oc r -4 (at « 30% of their virial 
radius), thus Hernquist profiles (Hernquist 1990) provide 
slightly better fits. 

• Mergers between halos in the cluster proximity occur 
with a frequency of about 5-10% since z=0.5. In the cluster 
environment mergers are rare; not a single merger occurs for 
halos whose orbits are contained within 1.6 Mpc = 80% -R200 
from the cluster center. 
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• The orbital distribution of substructure halos is close 
to that of an isotropic population of tracers in an isother- 
mal potential; the median value of apocentric to pericentric 
distances is 6:1, a value that does not vary with position 
within the cluster and is unbiased with respect to the orbits 
of the smooth particle background. Circular orbits are rare 
and about 20% of all our surviving halos within the cluster 
will pass within 200 kpc = 0.1i?2oo- 

• Most dark halos are tidally truncated to a value de- 
termined by the (average) density of the cluster at their 
pericentric positions. The approximation of isothermal halo 
mass distributions orbiting within a deeper isothermal po- 
tential works very well; i.e. rudal ~ r veri Ohaios / Ocius- 

• The mass attached to dark matter halos is approxi- 
mately 13% of the entire cluster mass and varies from 0% 
within ~ 200 kpc from the cluster center, to 20% at its virial 
radius. This latter value is roughly the expected value for the 
mass attached to halos above a circular velocity of 80 km s _1 . 
Correspondingly, the sizes of halos vary with cluster centric 
radius, an effect that may be observable using gravitational 
lensing of background galaxies. 

• Overmerging within the central regions of dense halos 
leads to a final distribution of substructure that is antibiased 
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